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ABSTRACT 

We present an analysis of data about Saturn's atmosphere from CassinVs 
Visual and Infrared Mapping Spectrometer (VIMS), focusing on the meteorology 
of the features seen in the 5-micron spectral window. We present VIMS mosaics 
and discuss the morphology and general characteristics of the features backlit 
by Saturn's thermal emission. We have also constructed a zonal wind profile 
from VIMS feature tracking observation sequences using an automated cloud 
feature tracker. Comparison with previously constructed profiles from Voyager 
and Cassini imaging data reveals broad similarities, suggesting minimal vertical 
shear of the zonal wind. However, areas of apparent wind shear are present in the 
VIMS zonal wind profile at jet stream cores. In particular, our analysis shows that 
the equatorial jet reaches speeds exceeding 450 m s -1 , similar to speeds obtained 
during the Voyager era. This suggests that recent inferences of relatively slower 
jet speeds of ~275-375 m s -1 are confined to the upper troposphere and that the 
deep (>1 bar) jet has not experienced a significant slowdown. Our measurements 
of the numerous dark, spotted features seen in the VIMS mosaics reveals that 
most of these features have diameters less than 1000 km and reside in confined 
zonal bands between jet stream cores. We propose that these spot features are 
vortices and that VIMS and ISS are sensing the same vortices at two different 
pressure levels. The local structure at the zonal jet streams remains complex, as 
VIMS may be sensing cloud features that are deeper than the NH 3 cloud deck. 

Subject headings: Saturn, Atmospheres — Atmospheres, Dynamics — Atmo- 
spheres, Structure 



-2 - 



Introduction 



Saturn's atmosphere has a reputation for being bland and unremarkable, especially when 
compared against Jupiter's colorful bands and vortices. A strongly scattering haze confined 
to the higher altitudes of the atmosphere is mostly responsible for obscuring the activity 
occurring at depth. Recent imaging results from Cassini have revealed this activity: dozens 



of compact, dark spots (likely anticylconic vortices) populate Saturn's mid-latitudes (Porco 



et al. 2005), and vigorous storms with inferred lightning discharges imply energetic updrafts 



in tandem with strong precipitation (Dyudina et al. 2007). Vasavada et al. (2006) describe 



tilted cloud streaks, abundant vortices, patchy bright clouds, and dark bands coincident with 
the location of zonal jet streams throughout Saturn's southern hemisphere. Their observed 
vortices had lifetimes on the timescale of several to hundreds of days, demonstrating the 
activity of the atmosphere. Though Cassini'' 's Imaging Science Subsystem (ISS) was not 



yet able to fully observe the northern hemisphere, Vasavada et al. (2006) noted similarities 



between both hemispheres based on previous studies from Voyager Smith et al. (1981); 



Smith et al. (1982); Sromovsky et al. (1983) 



Saturn's zonal wind profile, especially at its equatorial region, has been the subject of 
much discussion. Voyager imaging data indicated that the broad equatorial jet flowed as fast 
as 470 m s _1 at its peak (Sanchez-Lavega et al. 2000). However, Hubble Space Telescope 



(HST) observations of the jet from 1996-2002 showed evidence for a dramatic slowdown 



in the velocity of the jet by almost a factor of two {Sanchez-Lavega et al. 2003) over a 
broad (~30°) span of latitudes. Explanations for the slowdown include a true decrease in 
the speed of the jet or that HST observed higher altitude clouds that would flow at slower 
speeds ( Flasar et a/.|2005 ). Sayanagi and Showman (2007) performed numerical experiments 
that suggested that a combination of both is necessary to account for the apparent slowdown. 



Recent results from the Cassini orbiter made by Porco et al. (2005) have bridged the Voyager 
and HST measurements, as they supply a measurement of winds near 375 m s _1 at the 
equator using a near-infrared continuum band that samples deeper in the atmosphere than 
the HST measurements. 

With the arrival of Cassini in 2004, planetary scientists gained the capability of using 
an imaging spectrometer to assess the characteristics of the Saturnian atmosphere. The 
Visual and Infrared Mapping Spectrometer (VIMS) instrument is a 64x64 pixel spectrograph 
capable of near-simultaneous observations of a target from 0.3-5 /im VIMS is a hybrid 



instrument combining a slit-scanning visible spectrometer (Capaccioni et al. 1998) with a 



spot-scanning near-infrared spectrometer (Brown et al. 2004). For this paper, we focus 



primarily on data in the near-infrared, specifically at 5 /xm. We utilize this spectral window 
because of its potential to probe deeper layers of the atmosphere. In 5 /im, Saturn's thermal 
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emission can be directly observed, and Saturn's atmospheric cloud features are seen "backlit" 
against the infrared emission, indicating that these features are located anywhere in between 



the top of the atmosphere and the 5 fxm emission level. Bjoraker et al. (2007) estimates the 



5 jum emission layer at around 5 bars, which is below the estimated altitude (1-2 bars) for 
the cloud tops seen in ISS images observing reflected solar radiation. 



Baines et al. (2005b) and Baines et al. (2006) conducted an initial study of results from 
the first year of data from the Cassini VIMS instrument. Although their initial reports 
focused mostly on constraining atmospheric compositions, the reports made the important 
initial discovery of meteorological features seen in Saturn's atmosphere in the 5 fim spec- 



tral region. Momary et al. (2006) characterize these features, some of which have peculiar 



descriptions: certain features in the northern high latitudes are "donut-like" and are presum- 
ably vortices with a relatively cloud-less inner core. Furthermore, a train of cloud clearings 
arranged like a "string of pearls" inhabits the northern mid-latitudes. They determine that 
these clearings are quasi-evenly spaced with a near 4:1 flux contrast ratio and estimate their 



altitude to be near the 2.5 bar level. On a global scale, Momary et al. (2007) report a pro- 



nounced asymmetry in the hemispherical flux emitted by Saturn. The northern hemisphere 
emits nearly twice the flux of the southern hemisphere, possibly indicative of strong seasonal 
changes in the formation and dissipation of aerosols in the upper troposphere (~300 mbar) 
of Saturn. 

For this paper, we have undertaken an analysis of VIMS data of the Saturnian atmo- 
sphere from Cassini up to 31 March 2007. We present descriptions of the basic morphology 
of the features seen in the VIMS imagery, and provide statistics on the latitudinal and size 
distribution of these features. We also perform a comparative analysis of the features seen 
in VIMS and ISS mosaics. We have adapted an automated cloud feature tracking algorithm 
to construct the zonal wind profile of Saturn's atmosphere from VIMS data. We compare 
this profile to previously derived profiles to estimate the vertical wind shear. 

All latitudes in this paper, unless explicitly stated otherwise, are planetographic. 



2. Data and Methodology 



2.1. Data Sets and Common Procedures 



We obtained VIMS datacubes from the NASA PDS Imaging Node, and restricted our 
analysis to data sets of Saturn's atmosphere that either (i) attempted feature tracking by 
repeatedly observing the same absolute frame relative to System III longitude or (ii) moni- 
tored the overall state of Saturn's atmosphere via full-disk mosaics over a Saturnian day. A 
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table listing all of the data that we analyzed is included as Table [TJ We perform the standard 
VIMS data processing pipeline (e.g. removing cosmic ray hits, dividing by a pre-launch flat 



field, dividing by a solar spectrum, etc.) using software provided by the VIMS team. Barnes 



et al. (2007) discusses this pipeline further. We used the ISIS (Integrated Software for Im- 



agers and Spectrometers) software suite (Gaddis et al. 1997) for further data reduction and 



processing, which we discuss in the following paragraphs. 

All of the VIMS datacubes in this study were prepared for analysis by first co-adding 
the five longest wavelength channels in the data cube (5.057-5.122 /mi) and omitting the 
remainder. With ISIS, we re-projected the image into a rectangular (simple cylindrical) 
projection, and used kernels provided by the Navigational and Ancillary Information Facility 
(NAIF) at the Jet Propulsion Laboratory for absolute image navigation and orientation. Our 
re-projected images that were used for both mosaic construction and feature tracking are 
oversampled at 0.1° pixel -1 by utilizing a bicubic interpolation scheme supplied with ISIS. 
(1° on Saturn is ~1000 km.) Typical VIMS images have a native resolution of about 0.25° 
pixel -1 on Saturn, though this resolution strongly depends on the observation geometry and 
range from the target. Typical spacecraft ranges from Saturn during observations were ~9 
x 10 5 km, though these ranges sometimes varied by as much as 50%. 

A primary source of uncertainty in this paper is errors in the pointing knowledge of the 
spacecraft. Though the kernels provided by NAIF are carefully crafted, some uncertainty 
in the absolute location of where the spacecraft is pointing and the geographical location of 
the target is inevitable. We were unable to perform any manual correction to the pointing 
knowledge of the images via limb-fitting or other methods during the processing pipeline. 
We estimate the uncertainty in the absolute location of the cloud features to be 1-2° (J. 
Barnes, pers. comm.). However, the relative uncertainty in location between pixels is much 
smaller. Thus, our measurements of lengths and velocities, which are more dependent on 
relative uncertainty, do not degrade much in quality because of the small, but non-zero, 
absolute uncertainty. 



2.2. Mosaics 



We created global-scale mosaics of Saturn's atmosphere using VIMS data in order to 
characterize the frequency and size distribution of observed cloud features. In order to en- 
hance the visible features and maximize the data return from the VIMS instrument, we 



implemented an unsharp mask procedure described in he Mouelic et al. (2007). This tech- 



nique simply subtracts a percentage of a low-pass-filtered version of the original image from 
the original data. The resulting mosaics offer enhanced contrast and improved visibility for 
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the backlit features. No special processing was performed to mitigate image seams. 

We constructed two near-complete mosaics of Saturn's southern hemisphere, as Cassini 
observed the southern hemisphere twice over a period of nearly two weeks in December 2006. 
During both observation sequences, the spacecraft acquired data by imaging the hemisphere 
in a 3x3 grid (with one corner absent) at regular intervals during a Saturnian rotation. 
The entire southern hemisphere is observed up to 80° S, though some gaps in the coverage 
are present. We also constructed one mosaic of the northern hemisphere up to 60°N from 
observations in September 2006. This mosaic has relatively lower spatial resolution from 
unfavorable orbital and observational geometry. Other observational sequences that covered 
the northern hemisphere with improved resolution during the Cassini mission in the publicly 
available datasets were insufficient in global coverage or otherwise had unfavorable viewing 
conditions. 

2.3. Cloud Feature Census 

For the mosaics presented in this paper, we performed a statistical assessment of the 
population and characteristics of the features seen in Saturn's atmosphere. We limit our 
analysis to spot-like cloud features that are dark, compact, and ellipsoidal. These features 
are ubiquitous in the mid- and higher latitudes of the atmosphere. These features may 
be vortices, though we cannot explicitly confirm this because we currently are incapable 
of measuring local wind velocities within these structures. Furthermore, there are many 
complex and ambiguous structures where spots appear to be merging, shearing apart, or 
otherwise deforming. We elect to exclude these structures in our analysis in order to establish 
a lower bound in cloud feature population. To assess the statistics of the cloud features seen 
in 5 /im, we measured the central locations of these features to determine their latitudinal 
distribution. To measure lengths, we used the measuring tool in the ISIS software package 
to measure both the longitudinal and latitudinal span of these spots. This software tool is 
capable of precise, sub-pixel measurements, but we must be careful with their interpretation 
as we are examining over-sampled and interpolated mosaics. We estimate the uncertainty 
in the measurement of feature lengths to be at the level of 1 native resolution pixel (~250 
km). 
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2.4. Automated Feature Tracking 



We adapted our automated feature tracker ( Choi et aX]|2007 ) to measure the zonal wind 
profile of Saturn's atmosphere from VIMS cubes. Here, cloud features are assumed to act as 
passive tracers of the ambient zonal winds. However, instead of using hand and eye to track 
features in images separated by a discrete time, we use software to extract portions of an 
image and compare them to nearby portions of a later image. Our software then determines 
which portions have the highest correlation and measures their offset, which is then used to 
determine wind velocity. This technique has been used to measure winds of Jovian vortices 



Vasavada et al. (1998); Read et al. (2005); Choi et al. (2007)], eddy velocity fluxes on Jupiter 



and Saturn Salyk et al. (2006); Del Genio et al. (2007)], and Venusian atmospheric velocities 



{Toigo et al. 1994) 



We examined five particular data sets (Table [T]) out of more than twenty observation 
sequences that were designed for feature tracking. Unfortunately, we could not integrate 
the remaining data into our analysis due to poor observation geometry, map reprojection 
issues, poor image resolution, or a combination of all three. Our feature tracker normally 
attempts to measure velocities both in the zonal and meridional directions. However, we find 
no evidence of synoptic or regional scale meridional motion. Furthermore, the low spatial 
resolution of the VIMS data is insufficient for tracking local meridional motion. Thus, we 
modified our feature tracker to only measure zonal velocities. Subtle meridional shifts are 
apparent when manually blinking some images separated in time but are likely caused by 
uncertainties in camera pointing and image navigation. In addition, we adapted our feature 
tracker for the VIMS data by changing the size and shape of the correlation boxes that are 
extracted from the images. Small squares (the mode that our software typically uses) for 
correlation portions were unfit for the VIMS data as typical tracked cloud features were 
relatively large (a result of the low image resolution), and chances of returning a spurious 
result by the algorithm were high. We instead used correlation boxes that spanned a narrow 
range of latitudes (~1°) but the full range of longitudes available in the image. The technique 
returns a single zonal velocity for each latitude. The precise latitudinal location of each 



velocity vector is set at the mid-point latitude of each rectangular correlation box. Limaye 



(1986) used a similar technique to construct one of the first zonal wind profiles of Jupiter. 



Typical image separation times for VIMS cube frames that we analyzed range from 10 to 30 
minutes. 

The 1° latitude correlation box leads to a ~0.5° uncertainty for the vector's location, 
since the vector could result from feature(s) anywhere within the box. We can estimate the 
effect that the locational uncertainty has on the velocity measurement by multiplying it by 
the magnitude of latitudinal wind shear present. We estimate that the typical uncertainty 
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is ~10 m s -1 , though it could be up to 30-40 m s -1 in a high latitudinal shear environment. 
However, this is an upper bound to the uncertainty, as it is unlikely that the tracked features 
in a correlation box are all located at one edge of the box. Overall, the uncertainty in the 
zonal wind measurements is insufficient to significantly affect our conclusions. 



Results 



3.1. Atmospheric Features 

Much like the southern hemisphere, near-infrared continuum filter mosaics composed 



from Cassini's Imaging Science Subsystem (ISS) (Porco et al. 2004) analyzed by Vasavada 
et al. (2006), the VIMS mosaics (Figures flT[3j) exhibit a gradation with latitude in cloud 



feature clarity: the diffuse and muddled features in the lower latitudes contrast sharply with 
the more well-defined, distinct features seen in higher latitudes. The relative lack of well- 
defined, high-contrast features in the equatorial region supports the presence of an optically 
thick haze in the upper atmosphere capable of attenuating the thermal emission. However, 
occasional dark, wispy cloud features are visible in the equatorial region (Figure [4]). These 
are likely clouds at deeper altitudes that perhaps erupt periodically and are more effective at 
blocking the thermal emission. Equatorial 5 /mi hotspots similar to those encircling Jupiter 
(Ortiz et al. (1998) and Arregi et al. ( 2006[ )) have not been observed at low-latitudes, a 



result noted previously by Yanamandra- Fisher et al. (2001) through ground-based telescopic 
observations. 

The overall cloud morphologies in the VIMS and ISS mosaics are also broadly similar. 
The tropical latitudes in both hemispheres are dominated by numerous tilted, alternating 
dark and light linear features that strongly resemble the "tilted cloud streaks" noted by 
Vasavada et al. ( 2006[ ) in their study. An especially strong wind shear environment (10 to 
20 m s" 1 per degree of latitude) is present at these latitudes. Interspersed between the tilted 
features are zonally aligned light and dark stripes located at approximately 16°S, 22°S, 28°S, 
and 15°N. These bands span a wide range of longitudes but appear interrupted at some 
places and are not as distinct as the bands seen at higher latitudes. The region between 
30°S-35°S is a transitional region containing both occasional spot features (including a few 
with spiral patterns surrounding them, suggesting circulation) and stripes (both light and 
dark) occasionally interrupted by ambient clouds. A pair of narrowly separated zonal jets at 
31°S and 33°S are perhaps responsible for the complex meteorological patterns. A similar 
transition region containing two light stripes with ragged cloud features interspersed in 
between them is present between 30°N and 35°N. However, no known zonal jet exists at 
this latitude in the northern hemisphere. 
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The mid-latitudes (poleward of 35° latitude) in both hemispheres of Saturn are no- 
table for the presence of numerous compact dark spot features, indicating regions of high 
cloud opacity from their efficiency at blocking Saturn's thermal radiation. The majority 
of these spots are singular and isolated, though many instances of complex structures (ap- 
parent mergers, possible interacting structures, and amorphous shapes) are present inter- 



spersed among the singular spots. As discussed in Section |3.3[ these dark features are not 
distributed evenly across latitudes but instead concentrate in dense populations between dis- 
tinct, zonally-aligned stripes. These stripes correspond to locations of zonal jets (discussed 



further in Section 3.2), unlike the subtle bands seen in tropical latitudes, which are not 
associated with the jets. The image resolution and temporal coverage of the observations 
are insufficient for estimating the vorticity of these features, and we cannot definitively state 
how many of these spots are vortices. 

A salient feature of the northern hemisphere in the VIMS mosaics is the so-called "string 
of pearls" that exists in a narrow band just south of 40°N. These features of 
uncharacteristically high thermal emission, indicative of low cloud opacity. We catalogued 
22 "pearls" in the northern hemisphere mosaic, and examination of the central locations of 
these pearls reveals that they roughly divide into three groups. The first group of five and the 
second group of ten have similar spacings between pearls of about 3.7° longitude on average. 
However, the third group of seven are each spaced slightly wider, with a mean separation of 
approximately 4.8°. The shapes of the pearls are slightly elliptical, with a median aspect ratio 
(defined as east- west length divided by north-south length) of approximately 1.1. However, 
we note that there are some examples of pearls with shapes at both extremes, with aspect 
ratios ranging from 0.5 to 1.7. 

Two possible analogues currently exist for these features. Our preferred hypothesis 
is that the pearls are a manifestation of a von Karman vortex street, a repeating pattern 
of vortices where cyclones and anticyclones alternate with one another. The contrasting 
regions of cloud opacities could correspond to alternating regions of relative vorticity, with 
5 fim bright regions being cyclonic and dark regions being anticyclonic. The proximity 
of the pearls to 40°N, a quasi-westward flowing band wedged between strongly eastward 
flowing jets (see Figure [5]) supports the vortex street model. The White Ovals, a system 
of anticyclones thought to be part of a von Karman vortex street on Jupiter ( | Youssef and 



Marcus 2003), were located in a similar position in Jupiter's zonal wind profile (a westward 
flowing band in between two eastward jet streams). We again caution that we cannot directly 
measure the local wind flow of these features and thus cannot confirm that they are vortices. 
However, preliminary numerical simulations of Saturn's mid-latitude atmosphere support 
our hypothesis of a vortex street model for the pearls (K. Sayanagi, pers. comm.). An 
alternative possibility is that these areas of enhanced emission result from cloud modulation 
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caused by wave motion, similar to the 5 /im equatorial hotspots on Jupiter, though the 
"pearls" appear to be more numerous and smaller than their Jovian counterparts. Whereas 
equatorial hotspots on Jupiter range in length from 5,000 to 10,000 km (Orion et al. 1998), 
the average length of a pearl on Saturn is just over 1,000 km. Another contrast between 
the two features is the observation that there are typically ~10 Jovian 5 /im hotspots that 
globally span the latitude band between 5°N and 10°N (Arregi et al. 2006). In contrast, 
Saturnian pearls appear to be confined to a span of longitude approximately 100° in width, 
and are more numerous. The leading explanation for Jovian hotspots is that they are a 



manifestation of the downwelling branch of an equatorially trapped Rossby wave Allison 



(1990 


); 


Showman and Bowling 


( 


2000 


); 


Friedson 


( 


2(105) 



This physical model may also 
extend to the Saturnian pearls, though the wave would presumably not be equatorially 
trapped. A hybrid model with a Rossby wave acting to trap vortices and preventing them 
from merging is also possible. Youssef and Marcus ( 2003[ ) have suggested such a model in 
the case of Jupiter's White Ovals throughout much of the 20th century. 



3.2. Zonal Wind Profile 



Figure [5] is our constructed zonal wind profile. For comparison, we have also included 



the zonal wind profiles constructed from Voyager (Sanchez-Lav ega et al. 2000) and Cassini 



( Vasavada et al. 2006 ) imaging data. Generally, our zonal wind profile is in excellent agree- 



ment with the shape of previous zonal wind profiles from Voyager and Cassini and with 



previous preliminary reports about the VIMS dataset (Baines et al. 2005a), indicating that 



the overall jet structure has remained broadly consistent. An additional work by Baines 
et al. (2005b) states a measurement of 390 ± 50 m s _1 for a feature at 8°S, which is in 



agreement with our results. The broad equatorial jet is clearly evident, and its latitudinal 
shear is also in reasonably good agreement with previous zonal wind profiles. Our analysis 
of the VIMS images provides evidence for rapid equatorial flow at the jet core: clouds just 
north of the equator flow at ~500 m s -1 , and perhaps even faster. This suggests that the 
near factor-of-two reduction in speed relative to Voyager measurements for the equatorial jet 
( Sanchez-Lav ega et a/.|2003 ) and the more modest reduction measured by Porco et al. (2005) 
could perhaps be confined to the upper troposphere and that VIMS may be sensing deeper 
levels of the jet at these latitudes. We express caution with these particular measurements of 
rapid flow at the jet core (i.e. the cluster of points with u > 500 m s _1 at ~1-5°N and 7°S), 
as the equatorial region typically lacks discrete, high-contrast cloud features that are easily 
trackable by our technique. However, we have performed manual cloud tracking of a small 
cloud complex embedded in the flow (seen earlier as Figure [1]) that support the measurements 
above 500 m s -1 . We note that it is possible that the highest-speed, tracked cloud features 
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are a local disturbance and not characteristic of the zonally-averaged flow. Regardless, it 
seems clear that the equatorial jet speed exceeds ~400-450 m s _1 at the pressure sensed by 
VIMS 5 fim images. 

A few higher-latitude jets are present in the VIMS data, though unfortunately we cannot 
derive results south of 60° S using our current techniques on the publicly available data set. 
There is remarkably good agreement among all three zonal wind profiles for the precise 
latitudinal location of the jets. However, differences exist in the peak speeds measured at 
the jet cores, with the exception of the jet near 45°N, where our derived speeds from VIMS is 
nearly identical to that from Voyager. Furthermore, there is no universal increase or decrease 
of the peak velocities across all jets between the VIMS profile and past profiles. The jet cores 
at 75°N and 50°S appear to be flowing faster than in the previous profiles. The measurement 
at 75°N is attributed to a cluster of clouds located north of a large dark spot feature, and 
these clouds appear to be embedded in a relatively wide westward-flowing jet. At 50°S, 
a cloud feature that is embedded in a dark stripe (presumably corresponding to a narrow 
jet stream) is the likely source for our measurement. Though manual inspection of these 
features support our automated measurements, it is important to note that it is unclear how 
closely the motion of these features represent the ambient zonally-averaged zonal flow. The 
measurement at 60° S also supports the presence of apparent wind shear compared with past 
profiles, though we feel this measurement is less convincing from lack of repeat observations. 
One of the starkest differences between the VIMS zonal wind profile and previous profiles is 
the jet at 65°N, which appears to be flowing at nearly half of the speed of the jet at that 
latitude as seen by Voyager. Upon further manual inspection of the feature tracking images 
at that latitude, the evidence for wind shear at this location is inconclusive because of the 
low image resolution. 



3.3. Spot Feature Statistics 



Our database of spot features consists of 578 and 572 entries for the two southern 
hemisphere mosaics, and 297 entries for the northern hemisphere. For the purposes of com- 
parison, we have also analyzed the southern-hemisphere ISS mosaic composed and examined 
by Vasavada et al. ( 2006[ ), their Figure 1. This mosaic is composed of images observed using 
one of the near-infrared (750 nm) continuum filters. We only analyze the compact, patchy, 



light spots (clouds) in the ISS mosaic, as Vasavada et al. (2006) have already reported their 



analysis on the darker vortices. (It is important to note here that features that appear dark 
in VIMS 5 fim images are simply more effective at blocking the thermal emission, whereas 
features that appear dark in ISS images are inherently lower in albedo.) We catalogued 687 
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light spots in the ISS mosaic, and used the same ISIS analysis tools throughout all four 
mosaics. 

Figure[6]is a histogram showing the number of spot features as a function of their latitude 
for all four of the mosaics that we have analyzed. We sorted the features into latitude bins 
1° in size, centered at each integer latitude. The histograms for the southern hemisphere 
VIMS mosaics (Fig. [6^, and |6}d) demonstrate that the features are primarily located in 
three bands that are demarcated by zonal jet streams. The latitudinal distribution of these 
features is remarkably similar in the two mosaics, with minor differences between the two 
mosaics attributable to irregularities in observational coverage and inherent variability in 
the populations and lifetimes of these features. A fourth band containing spot features is 
suggested at higher latitudes (> 75°S). However, because of the poor image resolution and 
the cutoff in the image beyond 80° S, we are probably under- representing the population of 
features in this latitude band. 

Figure [6^ is the number distribution of the light, high-albedo spots in the ISS mosaic 
as a function of latitude. For comparison, we have included the population of the dark, 



low-albedo spot features counted by Vasavada et al. (2006) as black segments to supplement 



our results. The features are clearly divided into four latitudinal zones largely similar to 
the distribution seen in our southern hemisphere VIMS mosaics, with some differences in 
the number of features seen in each latitude zone between the ISS and VIMS mosaics. 
South of 50°S, we detect more features in the ISS mosaics than in the VIMS mosaics, with 
substantially more in the southern polar regions. This is likely a result of improved image 
resolution and facilitated feature identification at the higher latitudes in the ISS images 
(rather than incomplete longitudinal coverage, as we will show in Figure [7]). 

Figure [6]i is the number distribution of dark spots in the northern hemisphere VIMS 
mosaics. The tropical latitudes contain occasional spots, with the numbers increasing grad- 
ually with latitude. A main group of spots is centered poleward of 40°N, with a secondary 
group of spots present spanning 52-59°N. A broad gap in the spot population at 45-50°N 
aligns with a strong prograde jet at that latitude band. This matches the basic structure 
seen in the southern hemisphere, where zonal jets mark latitudes with a deficiency of spots. 
However, this structure does not extend to the tropical latitudes, as gaps in the population 
distribution exist on the flank of the fast equatorial jet. 

The histograms shown in Figure [6] are simply the number of spots that occur in latitude 
bins 1 ° wide, as a function of latitude, in our mosaics. However, our longitudinal coverage 
varies from latitude to latitude; moreover, the circumference of any given latitude circle 
decreases with increasing latitude. It is thus useful to normalize the spot numbers to remove 
these effects. To do so, we calculate the number of spots N that occur in latitude bins 1 ° 
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wide per unit longitudinal distance. This is given by 



Z,1T J JX curv 

where N is the number of spots in a latitude bin, / is the fraction of mosaic coverage 
at that latitude (unitless), and R curv is the zonal radius of curvature. Figure [7] shows this 
quantity for all four mosaics. Because N compensates for the gaps in the VIMS mosaic 
coverages, any remaining differences between Figures [7^ and b with [7^ suggest that other 
factors account for the discrepancies in the spot counts. The larger normalized count in the 
ISS mosaic at high latitudes (poleward of 60°S) compared to the VIMS counts in these bands 
is likely caused by improved image resolution in these areas. The smaller normalized count 
in the ISS mosaic at mid-latitudes (37°S-48°S) is more puzzling. This could be indicative of 
increased activity or a more complex vertical structure confined to these latitudes. 

Figure [8] presents a number-size distribution for all four mosaics, including the spot 



feature population in Vasavada et al. (2006). The overwhelming majority of the spot features 
are below 1500 km in east-west diameter. However, a few features have diameters that exceed 
2000 km, with one feature in both southern hemisphere VIMS mosaics exceeding 5000 km 
diameter. The ISS southern hemisphere distribution (Fig. ^p) reveals a slight difference 
in the location of the peak size. Whereas the VIMS southern hemisphere size distribution 
peaks at about 600 km, the ISS distribution peaks at a lower size, about 400 km. This 
may be an artifact of enhanced resolution in the ISS mosaics and as a result, improved 
detectability of the smallest features. The peak in the size distribution at 400-600 km has 
implications for Saturn's Rossby radius of deformation, a length scale at which rotational 
(Coriolis) effects become comparable to buoyancy effects. One can expect vortices at a size 
near the deformation radius, as structures that develop from inverse energy cascade and 



geostrophic adjustment typically attain sizes near this length scale Cho and Polvani (1996); 



Polvani et al. (1994)]. If we assume that the majority of this population are small vortices, it 



suggests that the deformation radius is within a factor of two of ~500 km. This estimate is 



relatively small compared to estimates of the deformation radius on Jupiter {Showman (2007) 



reports 1000-3000 km based on his numerical simulations) and in the upper troposphere of 



Saturn (2000-8000 km from Read et al. (2007) via an analysis of the zonal wind, temperature, 



and molecular hydrogen para-fraction profile). 

We examine a feature's east-west diameter as a function of its latitude in Figure [9] 
and Table [2] Clearly, most of the spot features are below 1500 km in diameter, with few 
features above this threshold. We note a subtle decrease in the characteristic diameters of 
the features with increasing latitude in both hemispheres. This slight trend also suggests 
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that the deformation radius is a controlling factor in the formation of these spot features. 
The deformation radius is inversely correlated with the Coriolis parameter /; therefore, the 
radius will decrease with higher latitude. 



Figure [10] presents the aspect ratios for all of the features catalogued in our study. We 
define aspect ratio as the ratio of a feature's principal east-west diameter to its principal 
north-south diameter. Most of the features span a broad range of aspect ratios, as there is no 
obvious bias towards a particular shape for these features in the main population below 1500 
km diameter. However, we note that larger features (> 2000 km) exhibit a preference for an 
aspect ratio greater than unity (zonal elongation). Aspect ratios for visible spots have been 



measured before on Jupiter and Saturn. Mac Low and Ingersoll (1986) noted that aspect 



ratios tended to be greater than unity with increasing east-west diameters based on their 



analysis of spots on Jupiter from Voyager imagery. Li et al. (2004) demonstrate the same 



trend based on their analysis of Jovian spots from Cassini imagery. Both studies reveal that 
spots on Jupiter with diameters above 3000 km tend to be more zonally elongated. Saturn 



appears to also follow this general rule. Vasavada et al. (2006) also show that the dark 



low-albedo vortices on Saturn tend to have aspect ratios greater than unity, though we note 
that from both Saturn studies, there is a general deficiency in the feature population above 
3000 km compared with the population found in Jupiter's atmosphere. This property of the 
features in both Jupiter's and Saturn's atmospheres may be a manifestation of the /3-effect 
(the variation of the Coriolis force with latitude) acting to preferentially elongate structures 
that exceed a diameter of ~2000 km. 



Discussion 



Bjoraker et al. (2007) estimates that the 5 /im thermal radiation emanates from a 



pressure level of 5 bars. However, the depth at which the thermal radiation originates does 
not necessarily have to coincide with the depth of the imaged features in the VIMS mosaics. 



Preliminary radiative transfer analysis by Baines et al. (2005a) and Momary et al. (2006) 



indicate that the features observed by VIMS are near the 2-bar level. Utilizing information 
gleamed from the ISS and VIMS zonal wind profiles, assuming that they originate from 
different depths, yields constraints about the latitudinal temperature gradient of Saturn's 
atmosphere. Saturn's winds are in geostrophic balance, and vertical shear of geostrophically 



balanced winds is controlled by the thermal wind equation (Holton 2004) 



dv 9 
dp 



4^ x V p T > 
Pf 
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where v 5 is the geostrophic wind, p is pressure, R is the specific gas constant, / is the 
Coriolis parameter, and T is the temperature. Thus, with all other factors being equal, ver- 
tical shear in the geostrophic zonal winds is directly correlated with the latitudinal gradient 
of temperature. If we assume that the features tracked by ISS and VIMS are separated by 
a scale height, we can place a bound in the vertical wind shear using the uncertainty in the 
zonal wind profile. Using our uncertainty estimate of ±10 m s -1 , dT/dy has an upper bound 
of ~0.5 K/1000 km. When extrapolated over 5,000-10,000 km, our zonal wind profile rules 
out latitudinal temperature contrasts exceeding 2-5 K throughout most of the atmosphere. 

However, apparent wind shear is present at certain jet stream cores. Our zonal wind 
profile shows faster wind speeds at the equator and at the 48°S jet core compared with 



previous profiles. Recently, Sdnchez-Lavega et al. (2007) used methane and near-infrared 



filter images from Cassini ISS, and observed higher-altitude features flowing at a slower 
speed compared to features at the main cloud deck (~700 mbar), yielding a vertical wind 
shear du/dz ~ 40 m s _1 per atmospheric scale height H (~50 km at the equator). Garcia 



Melendo et al. (2008) has extended that work to show that vertical wind shear (estimated 
at 10 m s" 1 per atmospheric scale height H) is present at several southern hemisphere zonal 
jets. However, if we again assume that the features tracked by ISS and VIMS are separated 
by a scale height, the ~100 m s -1 difference at the jet cores implies a shear that is greater 



than the estimate from Sdnchez-Lavega et al. (2007) by a factor of 2-3. (Our wind shear is 
also well above the ~25 m s -1 per atmospheric scale height at the 500 mbar level estimated 



by Flasar et al. (2005).) Importantly, our VIMS cloud tracking, coupled with ISS cloud 



tracking results, allows us to constrain the latitudinal temperature gradients at pressures 
significantly deeper than that directly measurable from thermal-emission spectra with, for 



example, Cassini's composite infrared spectrometer (CIRS) (Flasar et al. 2004) instrument. 
Our measured vertical wind shear implies a local latitudinal temperature gradient at the 2 
bar level to be ~1 K/1000 km at the equator, and as high as ~5 K/1000 km at the 48°S jet. 
Our estimates are upper bounds, as shear would be lower if VIMS were observing features in 
a deeper (> 2 bar) cloud deck. Thus, our analysis implies that the latitudinal temperature 
gradients at specific latitudes may exceed gradients found at other latitudes at the same 
pressure. 

From our analysis of the cloud features in the VIMS and ISS mosaics, it appears plausible 
that VIMS is simply observing the same cloud features that have been observed by visible 



imagers from Cassini and Voyager. Figure 11 supports our hypothesis. On the left of Figure 
11 is a portion of an ISS mosaic analyzed by Vasavada et al. (2006), whereas the right portion 



of Figure 11 is an inverted portion of Figure [TJ one of the southern hemisphere VIMS mosaics. 
Note the general similarities between both mosaics in the alignment of the horizontal stripes, 
the bright, white cloud features at higher latitudes, and the possible long-lived vortex just 
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north of 60°S. The appearance of the stripes at higher latitudes also match well, indicating 
that the high-albedo, white bands (suggestive of a thick cloud band) observed in ISS images 
are dark in the 5 fim images (indicating high opacity), and dark bands are bright in the 5 
ixm. data set (indicating low opacity). 

Comparison of the statistics of the spotted features in the VIMS and ISS mosaics also 
suggests that both instruments are typically observing the same features. Similarities in the 
features' latitudinal distribution and numerical population are evidence that the majority of 
the VIMS dark spots are white cloud patches seen in the ISS images. These bright cloudy 
areas likely contain abundant aerosols that are efficient at blocking abyssal thermal radiation. 
Unfortunately, this creates a degeneracy for the VIMS dark spots as both the patchy white 
clouds and the dark, low-albedo vortices contribute to the VIMS dark spot population. One 
key difference based on Figure [8] is that features analyzed by Vasavada et al. (2006) that 



exceed 1000 km diameter seem to be largely missing in the VIMS mosaics, as only a small 
percentage of VIMS features have diameters above this threshold. However, we note that 



based on Figure 8 of Vasavada et al. (2006), most of the dark vortices in the ISS images 
with diameters below 1000 km are simply classified as "dark." In contrast, most of the ISS 
vortices with diameters above 1000 km are classified as "dark with bright margin," "bright- 
centered," or simply "bright." Thus, it is possible that these more complex vortex families 
have a different structure when observed in 5 /im and were overlooked in our analysis. 

Overall, we believe that the majority of the dark spotted features seen in the VIMS 
mosaics are the same features as the ubiquitous white and dark spots seen in ISS images. 
In one scenario, ISS and VIMS could simply be observing the same cloud layer; lateral 
modulation of that cloud deck would then produce lateral variations in reflected sunlight 
(hence producing the ISS features) as well as lateral variations in 5-micron emission to space 
(hence producing the VIMS 5 fim features). However, estimates of the altitudes of features 
observed in ISS are ~700 mbar, which is substantially less than that suggested for the depth 
of the VIMS features (~2 bars). Alternatively, we propose a scenario where ISS is observing 
the reflected cloud tops of an ammonia cloud feature aloft, and VIMS is observing the base 
of this cloud feature at depth. This scenario implies that the features are nearly a scale 
height (or more) in thickness. For this scenario to work, the upper-level cloud tops that 
cause reflection of short-wavelength sunlight (as observed by ISS) must have small-enough 
particle size so that they do not block the 5 /im radiation upwelling from deeper levels. 
One thus needs a vertical gradient in particle size, with smaller particles aloft and larger 
particles at the base. The cloud structure associated with a given feature (e.g., a vortex) 
could thus cause reflection of sunlight from pressures of < 1 bar while allowing emission of 5 
/im radiation from ~2 bars. In any scenario, however, we must remain careful in interpreting 
the structure of Saturn's atmosphere at the jet streams, as VIMS may be observing tracers 
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of motion housed in a deeper NH 4 SH, or perhaps H 2 cloud deck. 



Though we have not directly observed rotation within these spot features, we believe 
that these spots are vortices. Our analysis suggests that the deformation radius of Saturn's 
atmosphere controls both their inherent size and the variance of their size with latitude. 
Furthermore, Penny et al. (2008) suggest that the Rhines effect ( Rhinesfl975 ) is suppressed 
at the latitudes containing VIMS spot features, indicating that these latitudes can generate 
and support vortices. Our proposed hypothesis that the spots have finite thickness is also 
consistent with the dynamics of coherent vortices; the vertical extent of such a feature is 
typically f/N times their horizontal width, where / is the Coriolis parameter and iV is the 
Brunt- Vaisala frequency Charney (1971); Reinaud et al. (2003)]. For Saturnian estimates 
(f ~ 10~ 4 s" 1 and N ~ 0.001-0.005 s" 1 ), a vortex 1000 km across would then be 20-100 km 
or 0.5-3 scale heights thick. Moreover, numerical models of vortex evolution show success in 



matching vortex observations when the vortices are 1-3 scale heights thick (Morales- Juberias 



and Bowling 2005). 



5. Summary and Future Work 

We have performed an analysis of an extensive data set about the Saturnian atmosphere 
in the 5 /zm spectral window from the VIMS instrument on board Cassini. Our mosaics of 
both the northern and southern hemispheres reveal an extensive population of dark, compact 
spots that are areas of blocked thermal emission. A statistical analysis of these features and 
comparison with earlier visible-light mosaics of the southern hemisphere suggest that the 
majority of the ordinary spot features seen in the VIMS mosaics are white cloud patches, 
with a minority population of dark (in reflected sunlight) vortices. Automated cloud feature 
tracking reveals that the general structure of the zonal wind profile is largely similar to 
profiles constructed from Voyager and Cassini data. Differences in the speeds exist at some 
latitudes, indicative of longitudinal variability or observation of deeper features indicating 
vertical wind shear. 

Many opportunities are available for continued analysis of the rich VIMS dataset con- 
cerning Saturn's atmosphere. Further radiative transfer studies will help constrain the alti- 
tudes of the features observed by VIMS and refine our estimates of the vertical wind shear 
and temperature gradients in the atmosphere. Analysis of features at polar latitudes was 
beyond the scope of this study, but interesting features are located at these latitudes, includ- 
ing a polar vortex above the south pole and a hexagonal standing wave above the northern 



polar regions (Baines et al. 2007). Cassini' s extended mission will yield an opportunity to 



measure seasonal changes in the chemistry, meteorology, and structure of the atmosphere. 
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Additional observation sequences designed for feature tracking will allow us to supplement 
our zonal wind profile and measure the degree of vertical wind shear throughout all latitudes. 
Furthermore, the majority of the VIMS data that we have analyzed are observations of the 
night side of Saturn, preventing us from comparing simultaneous VIMS and ISS images for 
common features. Comparative studies of the ISS and VIMS data sets and future planned 
observations will test vertical cloud structure models. We look forward to continued analysis 
of the data returned by VIMS and the Cassini spacecraft as we continue to explore Saturn's 
unique atmosphere. 
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Data set 


Latitude 


Spacecraft Clock Time 


Feature Tracking 


15°S-60°S 


1524828885-1524829781 




29°S-42°S 


1524977581-1524980093 




3°S-18°S 


1534393899-1534403233 




4°N-36°N 


1534575912-1534578687 




2°S-80°N 


1546355125-1546358126 


N. Hemisphere Mosaic I 


0°-60°N 


1536703543-1536741990 


S. Hemisphere Mosaic I 


0°-80°S 


1543657108-1543694887 


S. Hemisphere Mosaic II 


5°S-80°S 


1544708794-1544739494 



Table 1: VIMS data sets used in this study. 
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Fig. 1.— VIMS Southern Hemisphere Mosaic I. A snapshot of Saturn's southern 
hemisphere as observed in 5 /im by VIMS. This mosaic is an assembly of images taken on 
1 December 2006. Each frame component of the image has been projected in a cylindrical 
(rectangular) projection before final construction of the mosaic. 
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Fig. 2. — VIMS Southern Hemisphere Mosaic II. Same as Figure[TJ but approximately 
two weeks later. 
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Fig. 3. — VIMS Northern Hemisphere Mosaic. Saturn's northern hemisphere as seen 
by VIMS in the 5-micron spectral window. This mosaic is an assembly of images taken in 
April 2006. 
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Fig. 4.— Rectangular projection of VIMS data cubes in the 5-micron spectral window 
observing the northern tropical latitudes. The time separation of these two images is 29m 
22s. 
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Fig. 5. — Zonal wind profile for Saturn's atmosphere constructed from an automated feature 
tracker analysis of VIMS data. Our VIMS results are shown as small plusses. For comparison, 
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Uncertainties for the VIMS measurements are ~10 m s 1 , and are dependent on the amount 
of latitudinal wind shear around the measurement. 
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Fig. 6. — Histogram of the number of spot features as a function of latitude (planetographic). 
a. VIMS Southern Hemisphere Mosaic I. b. VIMS Southern Hemisphere Mosaic II. c. ISS 
Southern Hemisphere Mosaic, captured during the first orbit of Cassini around Saturn in 
2004. The unfilled areas are light spots analyzed in this current study, whereas black bars 
are dark, low-albedo spot features analyzed by Vasavada et al. (2006). d. VIMS Northern 
Hemisphere Mosaic. 
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Fig. 7. — Histogram of the number of spot features per 10,000 km as a function of latitude 
(planetographic). a. VIMS Southern Hemisphere Mosaic I. b. VIMS Southern Hemisphere 
Mosaic II. c. ISS Southern Hemisphere Mosaic, captured during the first orbit of Cassini 
around Saturn in 2004. (Both light and dark spots are represented in plot c.) d. VIMS 
Northern Hemisphere Mosaic. 
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Fig. 8. — Histogram of the number of spot features as a function of their horizontal (east- 
west) span. For clarity, we note the number for small bars that are difficult to distinguish, 
a. VIMS Southern Hemisphere Mosaic I. b. VIMS Southern Hemisphere Mosaic II. c. 
ISS Southern Hemisphere Mosaic. The black bars show the contribution from the features 
analyzed by Vasavada et al. (2006). d. VIMS Northern Hemisphere Mosaic. 
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Fig. 9. — Scatter plot showing a spot's east-west span as a function of its latitude (planeto- 
graphic). a. VIMS Southern Hemisphere Mosaic I. b. VIMS Southern Hemisphere Mosaic 
II. c. ISS Southern Hemisphere Mosaic, d. VIMS Northern Hemisphere Mosaic. 
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VIMS S. Hem I 


VIMS S. Hem II 


ISS S. Hem 


Latitude 


N 


Avg. Diameter (km) 


N 


Avg. Diameter (km) 


N 


Avg. Diameter (km) 


< 50°S 


196 


734.9 ± 203.5 


188 


718.6 ± 188.3 


89 


598.6 ± 170.4 


50°S-62°S 


180 


634.0 ± 200.5 


211 


609.4 ± 182.3 


179 


462.0 ± 139.0 


62°S-75°S 


187 


549.8 ± 175.2 


150 


525.0 ± 166.0 


272 


414.2 ± 113.9 


> 75° S 


14 


455.3 ± 109.6 


19 


527.0 ± 198.3 


146 


377.8 ± 145.3 


All 


577 


636.7 ± 207.9 


568 


620.5 ± 196.0 


686 


442.8 ± 151.2 





VIMS N. Hem 


Latitude 


N 


Avg. Diameter (km) 


20°N-30°N 


24 


937.2 ± 212.5 


30°N-37°N 


33 


977.5 ± 242.8 


37°N-46°N 


136 


860.5 ± 223.0 


46°N-60°N 


90 


817.0 ± 223.8 


All 


283 


866.8 ± 229.5 



Table 2: Average east-west diameters for the spot features in the northern and southern 
hemisphere mosaics, divided by latitude groups. Features with east-west diameters greater 
than 1500 km have been excluded from this analysis. 
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Fig. 10. — Plot comparing a spot feature's east-west span versus its north-south span. A 
line representing an equal aspect ratio is drawn for reference in each plot. a. VIMS Southern 
Hemisphere Mosaic I. b. VIMS Southern Hemisphere Mosaic II. c. ISS Southern Hemisphere 
Mosaic, d. VIMS Northern Hemisphere Mosaic. 
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Fig. 11. — (left) Portion of the southern hemisphere ISS mosaic published as Figure 1 in 
Vasavada et al. (2006). The contrast has been increased to enhance the observed features, 
(right) Inverted portion of the southern hemisphere VIMS mosaic seen in Figure [TJ 



